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Abstract. Thanks to advances in detector technology 
and observing techniques, true Sunyaev-Zel'dovich (SZ) 
surveys will soon become a reality. This opens up a new 
window into the Universe, in many ways analogous to the 
X-ray band and inherently well-adapted to reaching high 
redshifts. I discuss the nature, abundance and redshift dis- 
tributions of objects detectable in ground-based searches 
with state-of-the-art technology. An advantage of the SZ 
approach is that the total SZ flux density depends only 
on the thermal energy of the intracluster gas and not on 
its spatial or temperature structure, in contrast to the 
X-ray luminosity. Because ground-based surveys will be 
characterized by arcminute angular resolution, they will 
resolve a large fraction of the cluster population. I quan- 
tify the resulting consequences for the cluster selection 
function; these include less efficient cluster detection com- 
pared to idealized point sources and corresponding steeper 
integrated source counts. This implies, contrary to expec- 
tations based on a point source approximation, that deep 
surveys are better than wide ones in terms of maximizing 
the number of detected objects. At a given flux density 
sensitivity and angular resolution, searches at millimeter 
wavelengths (bolometers) are more efficient than centime- 
ter searches (radio), due to the form of the SZ spectrum. 
Possible ground-based surveys could discover up to ~ 100 
clusters per square degree at a wavelength of 2 mm and 
~ 10/sq. deg. at 1 cm, modeling clusters as a simple self- 
similar population. 

Key words: cosmic microwave background - Cosmology: 
observations - Cosmology: theory - large-scale structure 
of the Universe - Galaxies: clusters: general 



1. Introduction 

Cosmologists have long appreciated the value of the Uni- 
verse's biggest objects, galaxy clusters. Besides being a 
collection of galaxies well suited for studies of galaxy for- 
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mation, studies focussed on the global properties of clus- 
ters provide information on the nature of dark matter; 
the relative proportions of hot gas, dark matter and stars; 
and on scenarios of structure formation, including con- 
straints on the universal density parameter, Qq. One ex- 
ample of the latter that comes to mind in anticipation of 
observations with the new generation of X-ray satellites, 
Chandra and XMM, is the use of the redshift evolution of 
the cluster abundance to constrain ilo (Oukbir & Blan- 
chard 1992, 1997; Bartlett 1997; Henry 1997; BahcaU & 
Fan 1998; Borgani et al. 1999; Eke et al. 1998; Viana & 
Liddle 1999a); another is the now classic cluster baryon 
fraction test (White et al. 1993). 

While clusters have been extensively studied in the 
optical and X-ray bands, observations based on weak 
gravitation lensing and the Sunyaev-Zel'dovich (SZ) effect 
(Sunyaev & Zel'dovich 1972) are just coming to fruition. 
In the case of SZ observations, important samples consist- 
ing of several tens of clusters pre-selected in the X-ray are 
beginning to permit cosmologists to capitalize on the po- 
tential of combined SZ/X-ray observations (Carlstrom et 
al. 1996, 1999). FuU maturity of the field wiU be heralded 
by the realization of purely SZ-based sky surveys. In what 
we might refer to as the "SZ-band" , one can then imagine 
performing cluster science analogous to what is now done 
in the X-ray, e.g., the construction of cluster counts, red- 
shift distributions, luminosity functions, etc., all viewed 
via the unique characteristics of the the SZ effect. For ex- 
ample, several authors have emphasized the advantages 
of the SZ effect, over similar X-ray based efforts, to con- 
strain Qq via the cluster redshift distribution, as well as to 
study cluster physics out to very large redshifts (provided 
the clusters are out there, the very question of flo itself) 
(Korolyov et al. 1986; Bond & Meyers 1991; Bartlett & 
Silk 1994; Markevitch et al. 1994; Barbosa et al. 1996; 
Eke et al. 1996; Colafrancesco et al. 1997; Holder et al. 
1999). Such pure SZ surveys will be performed: the Planck 
Surveyor will supply an almost full-sky catalog of several 
thousand clusters detected uniquely by their SZ signal; 
and advances in both detector technology and observing 
techniques now offer the exciting prospect of performing 



2 



J.G. Bartlett: Sunyaev-Zel'dovich Surveys: Analytic treatment of cluster detection 



purely SZ-based surveys from the ground, with both large 
format bolometer arrays and dedicated interferometers. 

I discuss in this paper some aspects of the science ac- 
cessible to pure SZ surveys by examining the nature of 
their cluster selection. Because of the close analogy with 
X~ray studies, it is useful for this purpose to compare 
and contrast SZ-based cluster searches to those based on 
X-ray observations. The redshift independence of the sur- 
face brightness of a cluster (of given properties) means 
that SZ cluster detection is inherently more efficient than 
X-ray detection at finding high redshift objects. Equally 
important is that although the SZ effect and X-rays both 
"see" the hot intracluster medium (ICM), they do so in 
significantly different ways. In particular, the well-known 
fact that the SZ effect scales as the gas pressure implies 
that the flux density, S^, is simply proportional to the 
total thermal energy of the gas. This makes modeling es- 
pecially simple, for this quantity depends only on the to- 
tal gas mass and the effectiveness of gas heating during 
collapse, in stark contrast to the X-ray emission that de- 
pends also on the density and temperature distribution 
of the gas. This simplicity is an advantage because any 
theoretical interpretation of survey results requires an ad- 
equately modeled relation between the observable and the 
theoretically relevant quantity of cluster mass. 

These remarks concern essentially the physics of the 
'emission' mechanism itself. Of equal relavence is the na- 
ture of object selection imposed by the eventual detection 
algorithm used to extract sources from a set of observa- 
tions (a map); and this in turn depends crucially, as for 
any survey, on the particular combination of sensitivity 
and angular resolution of the observations. The objects 
detected by Planck will not be the same as those selected 
by ground-based surveys, and the final catalogs should be 
viewed as complementary. Planck will produce a shallow 
(~ tens of mjy) large-area survey, while the ground-based 
instruments will perform deeper surveys (< 1 mJy) over 
smaller sky areas (several square degrees). Most clusters 
remain unresolved at the Planck resolution of ~ 5 — 10 
arcmins, and this characterizes the kinds of objects acces- 
sible to this survey, e.g., the counts and the redshift distri- 
butions. The higher angular resolution of future ground- 
based instruments (on the order of an arcmin) will resolve 
many clusters and impose different selection criteria that 
will define the counts and redshift distributions of the final 
catalog. 

This is a central issue of the present study were, moti- 
vated by the possibility of ground-based surveys, I exam- 
ine the detection of resolved clusters. While the detection 
of unresolved sources is principally dependent on observa- 
tional sensitivity, and the final selection is more or less one 
of apparent flux ~ S,y ^ O^ii, - the detection of resolved 
sources is a more complicated cuisine involving individu- 
ally the characteristic source size, 9c^ and surface bright- 
ness, V. The specific goal of the present work is to quantify 
in terms of observing parameters the abundance, masses 



and redshifts of clusters detectable by ground-based sur- 
veys, with the particular aims of understanding optimal 
object extraction and the accessible science. For example, 
one of the key questions facing any survey is one of ob- 
serving strategy: given a fixed, total amount of observing 
time, should one "go deep" , with long integrations on a 
few fields, or instead "go wide" , covering more fields to 
higher sensitivity. If one is out to maximize the number of 
detected objects, the answer depends on the slope of the 
counts. One gains by going deeper if the integrated counts 
are steeper than S^"^, assuming that noise diminishes as 
1/Vt; otherwise, a larger area yields more objects. 

The cluster selection criteria of a survey may be com- 
pactly summarized by a minimum detectable mass as 
a function of redshift - Mdot(^)- Together with a suit- 
able mass function (we shall use the formalism of Press 
and Schechter 1974), this quantity determines both the 
source counts and redshift distributions of the final source 
catalog. Thus, in very concrete terms, we must examine 
Mdct{z) and understand the influence of the observation- 
ally imposed restrictions on 9c and i,^. Given a set of ob- 
servations, i.e., a map^ one could imagine many different 
algorithms to extract astrophysical sources, and Mdet(2^) 
will depend upon this choice. There is in principle an op- 
timal method, one which preserves signal-to-noise over 
the entire range of source surface brightness and size. It is 
characterized by a decreasing surface brightness limit with 
object size - the greater number of object pixels permits 
lower surface brightness detections. This algorithm is dif- 
ficult to apply in practice, and more standard approaches 
search instead for a minimum number of connected pixels 
above a preset threshold, thereby establishing a fixed cut 
on source surface brightness. Detection signal-to-noise is 
no longer constant, rather increasing with 0c , and these 
methods loose large, and in-principle detectable, low sur- 
face brightness objects. All of this will be reflected in the 
resulting functions M^ct{z). 

Throughout the discussion, we will be guided by the 
characteristics of two potential types of ground-based in- 
struments: large format bolometer arrays, epitomized by 
BOLOCAM (Glenn et al. 19980), and interferometer ar- 
rays optimized for SZ observations, as suggested by Carl- 
strom et al. (1999)0. BOLOCAM is a 151-element bolome- 
ter array under construction at Caltech for operation in 
three bands - 2.1 mm, 1.38 mm (the null of the thermal 
SZ effect) and 850 ^m. At the Caltech Submillimeter Ob- 
servatory, it is expected that the array will be diffraction 
limited to ^ arcminute resolution, or better, and limited 
in sensitivity by atmospheric emission (rather than de- 
tector noise). With its 9-arcmin fleld-of-view, one could 
imagine surveying a square degree to sub-mJy sensitiv- 
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ity in these bands. Carlstrom et al. (1999) have recently 
expounded the virtues of interferometric techniques using 
telescope arrays specifically designed for SZ observations. 
They have proposed the construction of such an array, 
operating at a wavelength of 1 cm, and estimated that it 
would be capable, in the course of one year of dedicated 
observations, of covering ~ 10 square degrees to a limit- 
ing sensitivity of ~ 0.3 mJy at arcminute resolution. In 
summary, then, we are interested in considering SZ obser- 
vations at arcminute angular resolution and to sub-mJy 
sensitivity at both centimeter and millimeter wavelengths. 

The paper is organized as follows: a rapid review of the 
SZ effect is given in the next section, followed by a discus- 
sion of the unique aspects of SZ cluster detection. Section 
3 details the cluster population model employed, based on 
the Press-Schechter (Press & Schechter, 1974) mass func- 
tion and the isothermal /3-model. Since we shall focus on 
issues of cluster selection as imposed by survey parame- 
ters, the cluster model will be restricted to the simple ex- 
ample of a self-similar population. The next section (Sec- 
tion 4) introduces the principal figures (Figures 1,2 and 3) 
of the present work by consideration of unresolved clus- 
ter detection; this case will also be used as a benchmark 
against which to examine the effects of resolved detection. 
Section 5 then develops the principle themes of resolved 
SZ cluster detection, starting with consideration of the 
optimal, constant signal-to-noise method, and followed 
by detailed study of cluster detection based on the stan- 
dard algorithm. A final discussion (Section 6) then more 
closely examines the number of detections to be expected 
from ground-based surveys and gives a non-exhaustive 
list of some important issues still to be treated. Section 7 
concludes. 

Key results will be the Mdct(-z) curves presented in 
Figure 1, quantifying the nature of SZ detected clusters, 
and the conclusion that resolved source counts are lower 
and steeper than expectations based on simple unresolved 
source count calculations, Figure 2. To the point, the latter 
implies that surveys at arcminute resolution gain objects 
with an observing strategy of "going deep" . The cosmolog- 
ical density parameter is denoted by VIq = SnGp/dH^, the 
vacuum density parameter by Ao = A/3H^ and the Hub- 
ble constant by Ho = hlQQ km/s/Mpc; unless otherwise 
indicated, h = 1/2 and Ao = 0. 

2. The Particular Value of the SZ Effect 

We begin by establishing our notation in recalling the ba- 
sic formulas of the SZ effect. The change in surface bright- 
ness relative to the unperturbed cosmic microwave back- 
ground (CMB), caused by inverse Compton scattering in 
the hot ICM, is expressed as 



(1) 



where x = hpV /kT^ is a dimensionless frequency expressed 
in terms of the energy of the unperturbed CMB Planck 



spectrum at To = 2.725 K (Mather et al. |1999D 
spectral shape is embodied in the function 

SkTof 
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= (2.28 X lO'^mJy/arcmin^)/^ 
while the amplitude is given by the Compton y-parameter 
kT 
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an integral of the pressure along the line-of-sight at posi- 
tion 6 relative to the cluster center. Here, T is the temper- 
ature of the ICM (really, the electrons), nie is the electron 
rest mass, Ue the ICM electron density, and ctt is the 
Thompson cross section. Planck's constant is written in 
these expressions as hp, the speed of light in vacuum as c, 
and Boltzmann's constant as k. These formulae apply in 
the non-relativistic limit of low electron (and photon) en- 
ergies; relativistic extensions have recently been made by 
several authors (e.g., Rephaeli 1995; Stebbins 1997; Challi- 
nor & Lasenby 1998; Itoh et al 1998; Pointecouteau et al. 
1998; Sazonov & Sunyaev 1998). The spectral shape of 
the distortion is unique, becoming negative at wavelengths 
larger than ~ 1.4 mm (relative to "blank" sky) and pos- 
itive a shorter wavelengths. This offers a way of clearly 
separating the effect from other astrophysical emissions. 

All of the physics is in the Compton y-parameter, an 
apparently innocuous-looking expression. In fact, it holds 
the key to all of the pleasing aspects of the SZ mechanism. 
First of all, the conspicuous absence of an explicit red- 
shift dependence is the well-known result that the SZ sur- 
face brightness is redshift-independent, determined only 
by cluster properties. This should be contrasted to other 
emission mechanisms which all experience "cosmic dim- 
ming" [l (X (1-1- z)~*] due to the expansion of the Uni- 
verse. This is countered in the SZ effect by the increasing 
energy density towards higher z of the CMB, the source 
of photons for the effect. 

Another very important aspect of the SZ mechanism 
resides in the fact that its amplitude is proportional to 
the pressure, or thermal energy, of the ICM. This appears 
most clearly when we consider the total flux density from 
a cluster, found by integrating the surface brightness over 
the cluster face: 

.kT{M,z) 



SAx,M,z) = jAx)D-\z) / dV 



-ne{M, z)(Jt 



(4) 

The integral is over the entire virial volume of the cluster. 
In this expression, Da{z) is the angular-size distance in a 
Friedmann-Robertson- Walker metric - 

n i \-o u-^ \^°^ + - + ^oz ~ 1)" 



f72(H-z)2 



= cH-^D{z) (5) 
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where I introduce the dimensionless quantity D{z). We 
see clearly that the final result is simply proportional 
to the total thermal energy of the ICM, / dVnT. 
This is extremely important, because it means that the 
SZ flux density is insensitive (strictly speaking, completely 
so for the total flux density and for fixed thermal energy) 
to either the spatial distribution of the ICM or its tempera- 
ture structure, making modeling much simpler than in the 
case of X-ray emission. Consider that in X-ray modeling 
one prefers the X-ray temperature over luminosity as a 
more robust indicator of cluster mass, but even the tem- 
perature has some sensitivity to the gas distribution, be- 
cause it is all the same an emission weighted temperature 
that is actually observed. We would expect the tempera- 
ture appearing in the second line of Eq. (^, which is the 
true mean electron energy, to demonstrate an even bet- 
ter correlation with virial mass than the observed X-ray 
temperature. Simple scaling arguments lead one to believe 
that this correlation should be T ~ Tviriai ~ M^/^(l -I- z), 
from which we deduce 



(6) 



where /gas is the gas mass fraction contributed by the ICM 
to the total cluster mass. 

The SZ mechanism therefore conveniently reduces all 
the potential complexity of the ICM to just its total ther- 
mal energy, cx /gas < T >. This quantity may nevertheless 
be influenced by several factors. For example, the gas mass 
fraction in Eq. has carefully been written as a general 
function of both mass and redshift. In simulations this 
quantity is most often constant, the majority of gas being 
primordial and simply falling into the cluster at formation. 
One could imagine other possibilities (e.g., Bartlett & Silk 
1994; Colafrancesco & Vittorio 1994) that would lead to 
a more important dependence on either mass or redshift, 
although metallicity arguments seem to require that most 
of the gas be primordial, at least in the more massive sys- 
tems (Metzler & Evrard 1994; Elbaz et al. 1995). While it 
appears from numerical studies that shocking during clus- 
ter formation efliciently heats the ICM to 80%- 100% 
of the virial temperature (Metzler & Evrard 1994; Bryan 
& Norman 1998), additional sources of heating could in 
principle change the temperature of the gas relative to 
that of the potential, i.e., T ^ rviriai- Such heating may 
not always produce the most obvious effects - remember 
that it is the total thermal energy of the gas that counts, 
and understanding the change of this quantity with heat- 
ing in a gravitational potential requires careful modeling. 
Although models studied so far do not lead to a strong 
effect (Metzler & Evrard 1994), we shall at times be dis- 
cussing rather low mass systems, for which these effects 
are poorly understood theoretically and observationally. 
Finally, the exact form of the virial temperature-mass re- 
lation depends in part on the dark matter proflle of the col- 
lapsing proto-cluster; once again, numerical experiments 
seem to indicate that this does not change too much, i.e.. 



one finds a good T-M relation with rather small scatter 
(Evrard et al. 1996; Bryan & Norman 1998). Putting all 
of this together, a relation of the form (||) between the 
observable, S^, and cluster mass appears quite reasonable 
and rather robust; and in any case, the modeling uncer- 
tainties are always easier to understand than in the case 
of X-rays, due to the all important insensitivity of the SZ 
flux density to spatial/ temperature structure of the ICM. 

The conclusion is that the SZ flux density should be a 
very good halo mass detector, in principle sensitive to all 
halos with significant amounts of hot gas and over a large 
range of redshifts. All of these remarks concern to a large 
extent the total SZ flux density of Eq. (^, and therefore 
apply primarily to situations where the clusters are unre- 
solved. It is still true that, even when a cluster is resolved, 
the SZ signal is proportional to the total thermal energy of 
the gas, but now only of that portion contained within the 
column defined by the beam. After first outlining the clus- 
ter population model employed, we shall tackle in detail 
the additional complexities introduced by resolved cluster 
observations. 

3. Modeling the Source Population 

The central ingredient of a model for the cluster popula- 
tion and its evolution is the mass function, n(M, z), which 
gives the number density of collapsed, virialized objects 
as a function of mass and redshift. The exact form of this 
function depends on the statistical properties of the pri- 
mordial density fluctuations. For Inflationary-type scenar- 
ios, in which these fluctuations are Gaussian, a reasonable 
expression for the mass function appears to be the Press- 
Schechter formula (Press & Schechter 1974) 



n{M, z)dM 



2 < p> 
M 



i/(M, z) 



dlncr(M) 



dlnM 



M ^ ' 



The quantity (p) represents the comoving cosmic mass 
density and u{Ad,z) = 5c{z)/cr[M,z), with 5c equal to 
the critical linear over-density required for collapse and 
cr{M, z) the amplitude of the density perturbations on 
a mass scale M at redshift z. Numerical studies ascribe 
rather remarkable accuracy to the simple expression of 
Eq. (0) (Lacey & Cole 1994; Eke et al. 1996; Borgani et al. 
1999), and we shall adopt it in the following. More explic- 
itly, 5c{z,no,Xo) and cr(M, z) = CTo(M) x (i:'g(z)/i:'g(0)), 
with Dg{z, rio, Ao) being the linear growth factor. It is es- 
sentially through Dg that the dependence on cosmology 
(f2o, -^o) enters the mass function, with Qq being the more 
important of the two as the dependence on Aq is relatively 
weak (see, e.g., Bartlett 1997 for a detailed discussion). 
This dependence on Qq in the exponent means that the 
cluster abundance as a function of redshift is a very sensi- 
tive probe of the density parameter (e.g., Oukbir & Blan- 
chard 1992, 1997), and is the motivation for many efforts 
in all wavebands to flnd clusters at high redshifts. As em- 
phasized by several authors (Barbosa et al. 1996; Eke et 
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al. 1996; Colafrancesco et al. 1997; Bartlett et al. 1998; 
Holder & Carlstrom 1999; Mohr et al. 1999), the SZ ef- 
fect is particularly well positioned in this arena (see also 
below). 

It is clear that the important theoretical variables are 
cluster mass and redshift. Although redshift is directly 
measurable, the mass appearing in Eq. (j7|) must be trans- 
lated into an observational quantity suitable for the type 
of observations under consideration. As mentioned above, 
one of the pleasant features of the SZ effect is the simplic- 
ity of this relation. Using the simulations of Evrard et al. 
(1996) to normalize the T — M relation, we can quantita- 
tively express the total SZ flux density of a cluster (e.g., 
Eqs. I & I) as 



5, = (34mJy/i8/3)/,(a;)/gasf^y' 



Anl(^) 



1/3 



M. 



5/3 
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{l + z)D--\z) (8) 



where the mass Mi5 EE M/W^ Mq refers to the cluster 
virial mass and /gas is possibly a function of both mass 
and redshift (see also Barbosa et al. 1996, but note that 
the definition there of D{z) differs by a factor of 2) . Evrard 
(1997) finds /gas = 0.06 while Mohr et al. (1998) 

find marginal evidence for a decrease in lower mass sys- 
tems (see also Carlstrom et al. 1999 for recent work based 
on SZ images); there is little information on any possible 
evolution with redshift at present. Other quantities ap- 
pearing in this equation are the mean density contrast for 
virialization, Anl(z, f^o, Ao) (= 178 for fio = 1, Ao = 0), 
and the dimensionless functions fu and D{z) introduced 
in Eqs. (|) and (|). 

Observations for which clusters are unresolved mea- 
sure this total flux density, and therefore this is all that is 
needed in order to calculate the unresolved source counts, 
as we will do in the next section. For resolved sources, 
on the other hand, the detection criteria are more com- 
plicated. Contrary to the point source limit, the details 
of the cluster SZ profile now play an important role. I 
will employ a simple isothermal /3-model to describe this 
profile: 



(9) 



The exponent a — 0.5(3/3 — 1), where (3 is the ex- 
ponent of the three-dimensional ICM density profile: 
n cx (1 -|- r^/r^)^'^'^/^, being the physical core radius. 
Local X-ray observations indicate that (3 ~ 2/3, a value 
I adopt throughout for the calculations. In this case, 
a = 1/2, a rather significant value, as will be discussed 
shortly. This profile will be assumed to hold out to the 
virial radius, i?v, of the cluster. 

The /3-profile of Eq.(^ is empirically described by j/o, 
a sort of central surface brightness (actually, it is yoju that 
has units of surface brightness, but it is simpler to work 
with j/o), and 9^- In these terms, there is nothing specific 
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Fig. 1. a)Detection mass as a function of redshift for un- 
resolved (dashed lines), optimal resolved (solid lines) and 
standard resolved (dot-dashed lines) detection satisfying 
gdotfpix = 1-5 mJy at a wavelength of 2 mm. In the unre- 
solved case, this simply corresponds to the limiting total 
flux density. For optimal resolved detection, g^ct refers to 
Qopt in Eq. (|lj), while for standard resolved detection it 
refers to ^st of Eq. ( p^ . The pixel size has been taken to 
be 0fwhm/2, and for the standard routine a detection an- 
gle 6'dot = 1/26'fwhm has been assumed, as indicated. In all 
cases these parameters correspond to 3f7 detections (see 
text for more detail). The upper (red) curves in each case 
correspond to the open model with f2 = 0.3. 

to the SZ effect. The physics of the SZ effect appears only 
when we make the connection between these empirical pa- 
rameters and the theoretically interesting ones, namely, 
mass and redshift, via relations of the kind yo{M, z) and 
9c{M, z). As our principle goal in this work is to under- 
stand the selection effects of resolved SZ cluster detection, 
the model for cluster evolution will be kept simple: a con- 
stant gas mass fraction, /gas = 0.06 h~^-^ (Evrard 1997), 
over cluster mass and redshift, and a core radius scal- 
ing with the virial radius R^, i.e., Xv = Rv/fc — const. 
Unless otherwise specified, this constant will be 
given a value of 10. One deduces from simple scaling 
arguments that 

1/3 

i?v = (1.69/i-2/3 Mpc) m'/-^(1 -I- z)-^n-^^^ 
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where the normalization is taken from the spherical col- 
lapse model. This scaling relation is about as robust as the 
relation for cluster temperature; in fact, the two are es- 
sentially the same, since T ~ M/R^. Some dependence of 
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Fig. 1. b) Detection mass at 1 cm for the listed parame- 
ters. The curves are labeled as in the previous figure, but 
note the change in scale along the y-axis. 

the normalization on mass and redshift could appear if the 
density profile around a peak forming a cluster changed 
significantly with these two quantities. In the following, 
we shall ignore this possibility, which numerical simula- 
tions seem to indicate is a small efi^ect in any case. This 
then fixes the relation 



rc(M,z) =i?v(Af,z)/a;v 



(10) 



For the axially symmetric surface brightness of Eq. 
the integral defining the total SZ flux density may be writ- 
ten 



s^iM,z) = j^27r J deey{e) 



1 + xi - 1 



Using Eq. (g) for Sv{M,z) in this expression, we deduce 

'Anl(2)' 



yo{M,z) = (6.40 X 10-"/i')/gasr!o 



Mi5(l + z)3 



178 



v/1 + 4 - 1 



(11) 



Together with the /3-profile (Eq. |), Eqs. (^ and 
(|ll|) define our cluster evolution model. As mentioned, 
it is self-similar, and we see the expected scaling 
Tc M^/'^ + z)-^ and yo M{\ + zf. This is most 
probably an oversimplified description of the actual clus- 
ter population, but it nevertheless provides a 'standard' 
with which we may understand the nature of the selection 



Fig. 1. c)Detection mass at 1 cm and for fipix = 1 arcmin^ 
(^fwhm — 2 arcmins) . The curves are labeled as in the pre- 
vious figures. Relative to Figure lb, the lower resolution 
results in smaller detection masses (note again the change 
in ordinate scale). The unresolved detection curves are 
unaffected by the change in resolution. 



effects imposed by resolved cluster detection, and a bench- 
mark for comparing more detailed models. It is important 
in the following that one does not forget the model depen- 
dence of our results, which can be retraced to this point of 
the discussion. 



4. Unresolved Detections 

This section is dedicated to the simple case of unresolved 
SZ detection, which will be used as a reference in the fol- 
lowing discussion of resolved detection. It also offers an 
introduction to the main figures. Figures 1, 2 and 3, sum- 
marizing the essential results of the present work. They are 
constructed for two representative cosmologies: a critical 
model fio = 1, and an open model (Ao = 0) with flo = 0.3. 
For the counts and redshift distributions of Figures 2 and 
3, I have used a CDM-like power spectrum with "shape 
parameter" fixed at F = 0.25; both models are normalized 
to the present day abundance of X-ray clusters - cs = 0.6 
and (78 = 1-0 for the critical and open models, respectively 
(e.g., Blanchard et al. 1999; Borgani et al. 1999; Viana & 
Liddle 1999b). 

Observations for which most clusters are unresolved 
measure the total SZ flux density. One can then simply 
invert Eq. (^ to find the corresponding limiting detection 
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mass as a function of redshift, M^jJ'^(z, S,y): 
JC,(».S„) = (0J2xl0"/,-"-'Me)('A.y'° 

Z36/5(z)(l + 2)-3/5 (12) 

Integrating the mass function over redshift and over 
masses greater than this limit directly yields the source 
counts: 

— {>S)-rdz-^r dM—(Mz) (13) 

The corresponding redshift distribution is simply obtained 
as the integrand of the z-integral. 

Figure 1 compares the various detection masses as 
a function of redshift for observations at 2 mm, e.g., a 
bolometer array, and at 1 cm, representative of an inter- 
ferometer; in each case the upper (red) curve corresponds 
to the open model. For the moment, concentrate only on 
the the dashed lines, which give the result for unresolved 
detection, Eq. (|l2|) , at a flux density of Su = 1-5 mJy. 
These curves remain unchanged from Figure lb to Ic, 
both at 1 cm but differing in angular resolution, because 
resolution is irrelevant for point sources (ignoring source 
confusion issues). Observe that in all cases the detection 
mass decreases with redshift beyond z ~ 1. This remark- 
able behavior is directly attributable to the fact that the 
SZ surface brightness is independent of distance. As al- 
ready emphasized, the distance appearing in Eq. ( |l2| ) is 
the angular distance and not the luminosity distance, a 
factor of (1 -I- z)^ larger. At high z the redshift depen- 
dence therefore scales as ~ z~^^^, one power coming from 
the assumed redshift scaling of the virial temperature and 
the rest from the decrease in angular distance (focusing) 
as ^ 1/z. A self-similar cluster model, implicitly assumed 
in this context by the constancy of /gas, thus predicts 
that SZ observations are more sensitive to objects at large, 
rather than intermediate, redshifts. This overall behavior 
would not change even if we broke the self-similarity with 
a declining gas mass fraction with mass; such a depen- 
dence could only modify the rate of decrease with z. On 
the other hand, an explicit decrease in /gas with redshift 
stronger than (1 -I- z)^-^ would cause M^^^. to actually in- 
crease with redshift. It is perhaps not so surprising that 
at close range, small z, the detection mass also drops; this 
is simply due to the increasing angular size of the object 
creating an increase in total flux density (the source is 
assumed to always remain unresolved in this discussion). 

From the difference between Figures la and lb,c, we 
see that, at a given sensitivity, the 2 mm observations 
probe farther down in mass. This is nothing more than the 
spectral shape of the SZ effect, described by the function 
ji, : the biggest decrement occurs precisely near 2 mm (the 



maximum emission of the effect is around 750/im). The 
resolved detection mass limits, to be shortly discussed, 
depend also on the angular resolution. 

Source counts for the two cosmological scenarios are 
given in Figure 2. These have been calculated using Eq. 
( p^ and the appropriate detection mass. In order to shed 
some light on the importance of low mass objects to these 
results, the counts are presented in pairs, one curve for 
a low mass cut-off of 10^^ Mq and one for a cut-off of 
10^^ Mq. Note that the x-axis denotes the pixel noise, 
(Tpix, and not a limiting source flux density; in the present 
situation of unresolved detections, this just means that the 
corresponding limiting flux density is qdct x Cpix. 

The first thing to remark from Figure 2 is the large 
difference between the two cosmological models. The pres- 
ence of clusters at high redshift in a low-density model 
shows up in the integrated counts, as confirmed by the 
corresponding redshift distributions shown in Figure 3, 
where the huge difference in cluster abundance at large 
redshift is evident. It is for this reason that the redshift 
distribution of SZ sources is a potentially powerful tool for 
constraining JIq (Barbosa et al. 1996, Bartlett et al. 1998). 
This is of foremost importance and represents one of the 
primary motivating factors behind this type of survey. 

This situation of unresolved sources applies in practice 
to missions such as the Planck Surveyor, as discussed, for 
example, by Barbosa et al. (1996) and Aghanim et al. 
(1997). The higher angular resolution of possible ground- 
based surveys calls for examination of resolved source de- 
tection. 

5. Resolved Detections 

In this, the principle section of this paper, we treat in de- 
tail the issue of resolved SZ cluster detection. The context 
will be one of arcminute resolution (pixel size) and sub- 
mjy sensitivity, as targeted by the up-coming ground- 
based instruments. It is worth being very explicit about 
the nature of the observations: the simplest case to imag- 
ine corresponds to that of an image produced by a bolome- 
ter array, such as BOLOCAM. In this case each point on 
the image, a 'pixel', represents a sample point of the sky 
brightness, as transformed by the optics of the observing 
system. The optical response may be divided into that 
of the telescope-plus-atmosphere (defining the projection 
of the sky onto the focal plane) and the optics proper to 
the detector (which act on the focal-plane image). There 
is a difference between bolometer arrays and the familiar 
example of a CCD camera working in the visible. For the 
latter, atmospheric seeing and telescope optics project the 
sky onto the focal plane by convolving with a Gaussian, 
and the camera itself then convolves this focal-plane im- 
age with a square top-hat, one centered on each pixel. 
The difference with a bolometer array lies in the fact that 
the CCD camera defines sharp, well-defined pixel bound- 
aries, while a bolometer array, with its set of cones, con- 
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volves the focal-plane image with something closer to a 
Gaussian. This means that, unlike CCDs, the pixels of a 
bolometer array 'overlap' in the focal plane. This has lit- 
tle consequence for the ensuing discussion, but it is all the 
same worth keeping in mind. 

This picture is not completely accurate when it comes 
to interferometers. Such instruments actually directly 
sample the Fourier transform of the sky. The result r^y 
often be modeled by a real sky image convolved with'Q.n 
effective, synthesized beam, but this beam lacks sensit^- 
ity on large scales, i.e., large spatial wavelengths on^e 
sky (short baselines). Thus, the effective beam cannot •eot 
be precisely a Gaussian, and it is especially important to 
correctly model the loss of response on large scales for 
extended objects such as clusters. For the ensuing discus- 
sion, I adopt the bolometer picture, applying it at times 
rather indiscriminately to characterize ground-based ob- 
servations; a future work will consider the details specific 
to interferometric observations (see also the recent work 
of Holder et al. 1999). 

For a bolometer array, the response of the entire optical 
chain (atmosphere-teleseope-detector) is often adequately 
modeled as a bi-dimensional Gaussian (if one is lucky, 
a symmetric one!), and for proper sampling, respecting 
Shannon, the sample period must be 2 - 3 times smaller 
than the beam FWHM. We will characterize a survey by 
the pixel size and sensitivity per pixel of its images - fipix, 
a solid angle, and Upix, a flux density. Note that because 
the pixels 'overlap' in the focal plane, what precisely is 
meant by fipix is the square of the separation between 
sample points, 0pix; the concept is a bit more ambiguous 
than in the case of a CCD camera. Thus, proper sampling 
means that the pixel scale ^Ipix = ^pix < ^fwhm/^- K is 
also worth explicitly remarking that, in the following, I 
assume that the noise is uncorrelated (from pixel to pixel) 
and uniform over the image. 

Given, then, a map of the SZ sky, we would like to 
understand how to extract clusters and the nature of the 
selection imposed by our technique. In addition to the ob- 
servational parameters rJpix and Cpix, this will depend on 
the form of the extended emission of the sources, a compli- 
cation avoided in the case of unresolved cluster detection; 
this represents an important difference between the two 
situations. Employing the /3-model introduced previously, 
Eq. (^, we see that a cluster SZ profile may be described 
by a characteristic central surface brightness, yo, and an 
angular size, 9c (the core radius). When couched in terms 
of the purely empirical parameters of rJpix, Cpix, j/o and ^c, 
we have before us a rather classic and well-known prob- 
lem of Astronomy. The only difference with galaxies in the 
optical is the form of the source profile. All physics spe- 
cific to the SZ effect itself appears only in the relation of 
the empirical source descriptors - {yo,9c) ~ to the theo- 
retically meaningful ones of cluster mass, M, and redshift, 
z. 



Resolved SZ Source Counts (Bortlett 2000) 

0000.00 i — — — 




S (pixel noise in mJy) 

Fig. 2. a) Cluster integral source counts at 2 mm as a 
function of map pixel noise for the two cosmological mod- 
els introduced in the text. The angular resolution and 
sampling correspond to the situation of Figure la. Unre- 
solved, optimal resolved and standard resolved counts are 
shown, respectively, as the dashed, solid and dot-dashed 
lines; the upper (red) curve in each case corresponds to the 
open model with = 0.3. For unresolved detections, the 
limiting source flux density is simply gdct x (pixel noise) . 
The light dotted lines in the background indicate the crit- 
ical slope of —2. The fact that the resolved counts are 
steeper than this value implies that, down to low noise 
levels, deep integrations yield more objects than wide and 
shallow ones. 

The procedure in the following is then always the same: 
quantify the detection algorithm in terms of fipix and tXpix, 
and then translate this, via the isothermal /3-model, into 
a Mdct{z;flpix, <ypix)- I employ a notation where the im- 
minently interesting independent variables of a function 
appear before the ";" , and parameterizing ones afterward. 
Thus, as written, the detection mass is primarily a func- 
tion of redshift, parameterized by the survey properties 
f2pix and ijpix. This function teaches us about the kinds of 
objects we detect, and leads directly to the survey counts 
and the redshift distribution of our clusters, via Eq. (|l3|). 
These latter quantities are the key indicators of the science 
content of the survey. 

This procedure will be applied to two source extrac- 
tion methods in the following, and the results compared 
to those for an unresolved SZ survey. We will refer to the 
first as "optimal detection" , because it extracts sources in 
such a way as to preserve the signal-to-noise across the 
entire range of detectable surface brightness and source 
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Fig. 2. b) Same as the previous figure for an observation 
wavelength of 1 cm; curve types have the same meaning 
as before (the labels 'Optimal' and 'Standard' have been 
removed for clarity). The resolution and sampling corre- 
spond to the situation of Figure lb. 
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Fig. 2. c) Same as Figure 2b (A = 1 cm) but now for 
i^pix — 1 arcmin^, i.e., the situation of Figure Ic. The 
smaller detection masses at this lower resolution result in 
higher counts when compared to Figure 2b. Note that the 
unresolved counts are the same here as in Figure 2b. 



size. This is achieved by lowering the surface brightness 
limit for large sources, possible due to the greater num- 
ber of covered pixels. The second method, routinely used 
by such packages as SExtractor (Bertin & Arnouts 1996), 
searches for a minimum number of connected pixels above 
a preset threshold. The important difference with the first 
technique is the imposition of a fixed surface brightness 
limit, independent of source size. The signal-to-noise of 
the detections is no longer constant, but increases with 
source size. This technique may be considered sub-optimal 
in the sense that it loses in-principle detectable low sur- 
face brightness sources, a fact well appreciated in the case 
of optical galaxy surveys. 

5.1. Optimal case 

Optimal detection selects all sources with a flux density 

(assuming spatially uncorrelated and uniform noise) where 
N is the number of pixels covered by the cluster and 
Qopt represents a threshold, say Qopt ~ 3 — 5; in fact, 
Qopt = S/N, the signal-to-noise of the detection. Notice 
also that, as advertised, the limiting surface brightness de- 
creases with object size: < >~ S,j/N ~ Qoptc^pix/VN. 
One extracts in this way all objects detectable at a given 
S/N, and for this reason we may refer to the method as 
optimal. The number of object pixels is simply found as 
N = 7r0^;j./r2pix, where ^vir = ^v/-Dang is angular virial 
radius. This permits us to express the detection mass as 



(0.19 X IO^^/i^-^Mq 

2\ 3/8 



3/4 



arcmm 



mJy 

(/i^/gas) ^ ^ 



+ (14) 

As written, this criteria uses an aperture correspond- 
ing to the full angular size of the object - Si, is under- 
stood to be the total SZ flux density in Eq. (||). For re- 
solved sources, one would like to chose an aperture which 
optimizes the signal-to-noise ratio of the detection. In- 
terestingly, a 3D gas profile close to r~^, corresponding 
to a SZ surface brightness y oc 0~^, results in a con- 
stant signal-to-noise with aperture radius. A /3-model 
with n cx (1 + r'^ /rl)-^'^/'^ and /3 - 2/3 exhibits this be- 
havior at large radii, for example: y{9) ~ (1 -f 0^/0^)"^/^. 
In this case, the signal-to-noise of a SZ detection increases 
from the center of the cluster image out to the core radius, 
Tc, beyond which it turns over to a constant out to the 
virial radius. The situation is different for X-ray images, 
where the surface brightness falls off more rapidly, diving 
under the background at large radii. From this we con- 
clude that the simple criteria given above provides in fact 
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Fig. 3. a) Redshift distribution of the integrated counts 
for a flux density of gdctCpix — 1-5 mJy at 2 mm. The 
parameters are the same as those in Figure 2a, and the 
hne types have the same meaning. The upper (red) curve 
in each case corresponds to the open model with Qq — 0.3. 
Note the large difference between the two cosmological 
models apparent in all cases. 



an optimum SZ detection (at least as long as P remains 
close to 2/3, as appears to be the case locally). 

The detection mass Eq. (^4|) is displayed in Figure 1 as 
the solid lines. Compared to the hypothetical point source 
results, observations resolving clusters are less efHcient 
at detecting clusters, especially at intermediate redshifts. 
This is easy to understand as the effect of distributing a 
given flux density over N pixels, each adding a noise with 
variance of dpix, resulting in a total noise level over the 



object image of vNapi-^- A point source, in contrast, is 
only subject to the noise of one pixel, Cpix. Hence, high 
resolution at fixed sensitivity "resolves out" a certain frac- 
tion of objects. The consequences for the source counts are 
clear and will be discussed shortly. These curves retain 
the same asymptotic behavior as before, namely a greater 
sensitivity to low masses at high redshift. Despite the fact 
that the object covers a larger number of noisy pixels as 
z decreases, the optimal method is able to take proper 
advantage of the greater total flux density to detect low 
mass objects locally, just as in the case of unresolved point 
sources. We shall see that this does not follow for the stan- 
dard detection routine (the dot-dashed lines), due to its 
additional surface brightness constraint (discussed below) . 
By comparing Figures lb and Ic, which differ only in their 
angular resolution, we note that for a given sensitivity, 
lower resolution observations are the more effective. This 



Fig. 3. b) Redshift distribution of the integrated counts 
for a flux density of gdctCpix = 1-5 mJy at 1 cm, and for 
the same parameters as in Figure 2b. 



is traceable to the fact that the flux density of a source is 
dispersed over fewer (noisy) pixels than would be the case 
at a higher angular resolution. This indicates that low res- 
olution observations at a given wavelength and sensitivity 
are to be prefered, at least for detection purposes. There 
is, however, a limit set by eventual source confusion. 

We have just seen from Figure 1 that low surface 
brightness clusters are "resolved out" at high resolution. 
This leads to overall lower counts that are also much 
steeper than the equivalent for unresolved point sources. 
Generally speaking, the unresolved counts do not deviate 
too much from a Euclidean law, (x , on the other 

hand, the resolved counts can be much steeper. The ex- 
amples shown in Figure 2 are in fact steeper than 5*"^, in- 
dicated by the dotted lines, down to essentially the faintest 
flux levels attainable in immediately foreseeable observa- 
tions. This is critical for optimizing an observing strat- 
egy with a fixed amount telescope time, T. Consider the 
common situation in which the final map noise decreases 
with integration time as l/\/t; then, the solid angle cov- 
ered in time T, with individual field integrations of du- 
ration t, scales with sensitivity as ~ T/t ~ tTpj^. Hence, 
if the integrated source counts are steeper than cTpj^, one 
gains objects by "going deep" , integrating longer on each 
individual field, rather than "going wide" , with shorter 
integrations covering a larger total solid angle. The im- 
portant conclusion to draw from Figure 2 is then that the 
way to optimize the number of detected objects in a sur- 
vey with arcminute resolution is by "going deep" , down 
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Fig. 3. c) Redshift distribution of the integrated counts 
for a flux density of gdotCpix = 1-5 mJy at 1 cm, and for 
the same parameters as in Figure 2c; in particular, for 
^fwhm = 2 arcmins. 

to the point were the counts begin to flatten out. In our 
examples, this does not occur until the very lowest flux 
levels deemed at present reasonable. It should be empha- 
sized that this conclusion rests on the results calculated 
here in the context of a self-similar cluster population. It 
is all the same suggestive and important in the fact that 
it is contrary to the conclusion one would draw based on 
unresolved source count calculations. The opposite holds 
for surveys with low angular resolution where the majority 
of sources remain unresolved, such as the Planck Surveyor 
observations. 

Finally, as to be expected from the "loss" of objects at 
intermediate redshifts, the redshift distribution for opti- 
mally selected objects lies under the corresponding point 
source examples, and is somewhat flatter. All the same, 
the two cosmological models are easily distinguished with 
an enormous "leverage" at high z. 

5.2. Standard algorithms 

Standard detection routines typically identify sources as a 
minimum number of contiguous pixels all above a preset 
threshold, usually qst times the pixel noise (Tpix- This is not 
the same criteria as above, in the optimal case, because we 
have now established a fixed surface brightness threshold - 
feCpix/f^pix ~ independent of object size (or luminosity). 
Previously, we allowed ourselves to lower this threshold for 
larger sources, in order to pick-up low surface brightness 
objects while maintaining a constant signal-to-noise; for 



this reason, it was an optimum detection algorithm. Here, 
the surface brightness is instead a flxed, while the signal- 
to-noise increases with object size as S/N — qstVN- A 
further difference is that the surface brightness cut im- 
poses a minimum detectable mass at z = 0. We obviously 
expect this method to detect fewer objects than the opti- 
mal approach. 

Consider application of the standard algorithm to a SZ 
profile, empirically described in the /3-model by j/o and 9^. 
Although our final goal is to understand the selection on 
mass and redshift imposed by the detection criteria, it 
is quite useful, firstly, to gain insight into the workings 
of detection in terms of tJo and 9c. As mentioned, what 
is actually recorded at each sample point (pixel), say by 
a bolometer camera, is the sky signal integrated over the 
beam B, which we will take to be axially symmetric. Thus, 
for a pixel at position h (a unit vector on the sphere): 

Sf'%h) = / dVL'i^{n')B{h-n') 



For our calculations, we shall furthermore adopt a Gaus- 
sian beam, so that a cluster appears as a /3-profile smeared 
by a Gaussian of dispersion db = ^'fwhm/ V8 In 2, 

where 9 is the angle from the beam axis and 0fwhm is 
the beam full-width at half-maximum. Notice that we 
take the image to be in flux density units. By placing the 
coordinate origin at the cluster center, so that now h is 
simply marked by the angular distance 6 from the origin 
(small angle approximation), the beam-smeared profile of 
a cluster may be written as 

s:''\9)^y,el3,g[9/9c;a^,/ec] 

explicitly separating out a dimensionless function 



Q{r\p) 



dx xe p ' 



2xr cos ( 



(1 



■ 2xr cos ( 



parameterized only by the ratio p = a\^l9c. The Heavyside 
function, O, cuts off the integral beyond the virial radius. 

It is this smeared profile of a cluster that is subject 
to the detection criteria that a minimum number of con- 
nected pixels, A'niin, must lie above the threshold QstCpix- 
This amounts to demanding that the object image above 
a flux density of QstCpix cover a minimum solid angle of 



A^ in in 



Let ^dot be the angular size of a cluster above 



the detection threshold, which may be calculated as the 
root of the following equation: 



ST 



7dct 



) — 'ZstO'pix 



(15) 



We will say that a cluster is detected if 6'dot is large enough 
to cover A^min pixels. In the present analytic treatment, 
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Fig. 4. a) The behavior of 6'dct/^'c as a function of 9c/ (Jh 
and Uo in a critical-density universe (f2 = 1, /i = 1/2) and 
for \y\ = 9.48 x 10^^, e.g., A = 2 mm, ^fwhm ~ 1 arcmin 
and gdct.CTpix = 1-5 mJy. 




Fig. 4. b) Same as Fig. 4a, but for \y\ = 6.84 x 10"'', e.g., 
A = 1 cm, 6'fwhm = 1 arcmin and QdctCpix — 1.5 mJy (or 
A = 2 mm, 6'fwhm — 1 arcmin and gdct^pix = 11 niJy, or 
A = 2 mm, 6'fwhm = 2.69 arcmins and gdotCpix = 1-5 mJy). 



we will simply impose a lower limit to ^dot and ignore any 
complications arising from the discreteness of the image. 
Using the function Eq. ( p^ may be written in compact 
form as 

Q[e,c,iec;a^l6c] = ^ ^-{t) (16) 

Vo"c3v Vo V "c / 

introducing the parameter y = (5stO'pix)/(o'^Ji/) character- 
izing the experimental set-up. It is clear that the solution 
will be given as ^dct/^c, that it will be a function of <t\^/9c 
and 2/o, and that it will be parameterized by y, i.e., 

[ddct/Oc]{(Jh/dc,yo;y) 

To understand the role of y, study the result in the limit as 
— oo; this will be particularly important below, when 
we consider the non-zero detection mass at zero redshift 
imposed by the surface brightness cut. In this large-object 
limit, g{r;p) g{r ^ Q;p ^ Q) ^ 2ttp^ . We thus find 
that in order to be detected an object must have a central 
surface brightness 

v. > I (17) 

In other words, y indeed embodies the surface brightness 
cut. This will be used shortly. 

Figure 4 shows the solution over the (^c/cb, yo)^plane 
for two reasonable values of \y\. To understand this fig- 
ure, separate the plane into a region occupied by resolved 
sources - 9c/<Jh >> 1 - and the region of point sources, 
9cl(J\, « 1: 



— Resolved sources: It is clear that by increasing j/o at 
fixed 6c/(Jh{» 1), we see an ever increasing portion 
of the ICM. The cluster 'lights-up' until we see all of 
it, out to the virial radius (beyond which we assume 
that the gas has not been heated), and the solution 
flattens out at this point to the adopted value cCv = 
Rv/fc — 10; beyond this, there is no more cluster to 
be seen. Of course, in the other direction, the object 
is eventually lost as we decrease yo to the point where 
even the central parts of the cluster do not rise above 
the detection threshold. 

— Point-source limit {9c/<Jh 0): In this extreme, the 
source profile becomes that of the beam, normalized to 
the total source flux density. This latter quantity scales 
as yojuOc, so that as ^c/^b continues to decrease at 
fixed y (i.e., holding constant), the imprint of the 
object gradually sinks below the detection threshold 
and ^dct — > 0; this explains the cut-off at low 6c/ for 
a given central surface brightness. 

So far, nothing extraordinary, but rather the standard 
issues of extended object detection given a particular in- 
tensity profile. Modeling more specific to the SZ effect en- 
ters only when we apply the SZ-based relations between 
the central surface brightness and angular size - j/o and 
6c - and the theoretically meaningful quantities of clus- 
ter mass, M , and redshift, z. These relations allow us to 
translate a surface like that of Figure 4 into an equivalent 
surface over the (z, M)-plane, as shown for three different 
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Fig. 5. a) Source detection radius (arcmins) as a function 
of cluster mass M (units of 10^^ solar masses) and redshift, 
for 6'twhm = 1 arcmin and y — 9.48 x 10~^, e.g., A = 2 mm 
and (^stCpix — 1-5 mJy. Notice the non-zero detection mass 
at z = 0. 

cases in Figure 5 (all using our adopted self-similar clus- 
ter model). Note that these latter surfaces (Figure 5) are 
not uniquely parameterized by y, because the translation 
from the axes in Figure 4 to mass and redshift explicitly 
involves (Jh- Thus, there are now two governing parame- 
ters, which we can take to be y and ai^: 

Odd [z,M-,y,ah] 

To study Figure 5 in detail, recall the simple scal- 
ing relations j/o ^ nM ^ AfANL(z)(l -I- z)^ and 9c ~ 
M^/^A^l^^{z){l + z)-^/D^ng{z), valid if the core radius 
scales with virial radius (it may not, but this has been 
assumed in the construction of the figure). Notice that 
mass and redshift are mixed in a nontrivial way in the 
expressions for surface brightness and angular extent. In 
particular, a cluster of given mass becomes more centrally 
bright towards higher redshift, due to a higher gas den- 
sity (scaling the background), while its angular extent 
at first decreases rapidly, as 1/z at low redshift, and 
then approaches an asymptote, since (1 -I- z)Dang(2) — > 
2c{Ho^o)~^ towards large z. 

Consider firstly the low-redshift region of Figure 5b 
(the various effects are most clearly displayed in panel 
b), where the z-dependence is dominated by Uang; here, 
mass uniquely parameterizes the surface brightness, i.e., 
j/o, while z changes only 6c'. 

— At constant M, decreasing z increases the angular 
size of a cluster, so that, for resolved objects that are 




Fig. 5. b) Source detection radius for 0fwhm = 1 arcmin 
and y ~ 6.84 x 10~^, e.g., A = 1 cm and ^stCpix = 1-5 
mJy. Note the change of scale for the M-axis relative to 
the previous figure. 

not completely below the surface brightness detection 
threshold (qstCpix/^^pix), ^'dct rises as 1/z; this corre- 
sponds to the Odet/dc const behavior in Figure 4. For 
smaller objects, of low mass, the beam profile deter- 




Fig. 5. c) Source detection radius for 0fwhm — 2 arcmins 
and y = 6.84 x 10^"*, e.g., A = 1 cm and (/stCpix = 1.5 mJy. 
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mines the scaling of ^dot with z, and this corresponds 
to the point -source limit of Figure 4. 

— At fixed (low) redshift, low mass objects eventually 
fall below the surface brightness limit, and ^dct reaches 
zero; on the other hand, massive clusters are 'lit-up' 
out to their virial radius, at which point ^dct attains 
^vir = a^v^c, which grows as M^/'^. 

An important aspect of detection procedure, mentioned at 
the beginning of this section and now evident from Figure 
5, is the existence of a minimum detectable mass in the 
limit of zero redshift (on the M-axis). This characteris- 
tic mass is established by the surface brightness threshold 
- feCpix/i^pix, and its existence represents a fundamen- 
tal difference with the previous case of optimal detection, 
where arbitrarily low mass (low surface brightness) clus- 
ters where picked up if they were large enough, i.e., very 
close at z = 0. We already saw in Eq. (0) how y summa- 
rizes the surface brightness constraint. The low redshift 
detection mass limits seen in Figure 1 and here in Figure 
5 are indeed reproduced numerically from Eq. (p^, once 
Eq. (|ri|) is used to convert j/o into a mass at z = 0. 

At redshifts approaching unity and beyond, M and z 
are fully mixed in the expressions for j/o and 9c'- 

— Massive clusters well above the surface brightness 
threshold increase in surface brightness with redshift 
to the point where they are completely seen, all the 
way out to 6'vir; at even larger z, Odot reflects the grad- 
ual fall-off to the asymptote set by the angular-size 
distance. This explains the ridge running down the sur- 
face in Figure 5 around z = \. Less massive clusters, 
on the other hand, only reach the point of full illumi- 
nation at higher z, well into the asymptotic behavior 
of 0vir, and hence the ridge tends to be washed out at 
the low mass end. Finally, the central surface bright- 
ness of very low mass clusters falls below the detection 
threshold at ever larger redshifts, i.e., the boundary 
^det = moves outward in z as M is decreased. 

Compare now the three panels of Figure 5. We observe 
the greater sensitivity at 2 mm, due to the spectrum of 
the SZ effect, by the fact that a given cluster of mass M 
and z produces a smaller ^dct at 1 cm wavelength in panel 
b (notice the change in scale along the M-axis between 
panels a and b). The same remarque applies to the greater 
sensitivity, at a given noise level, of the lower resolution 
observations exemplified in panel c. These characteristics 
will be inherited by the detection mass curves, our next 
topic. 

5.2.1. Detection mass as a function of redshift 

Since O^ct increases monotonically with Af, the contours 
displayed on the top and bottom faces of Figure 5 rep- 
resent curves of minimum detectable mass, M|J;^(z; ^dct), 
each one for a different detection threshold defined by dif- 



ferent values of 0dct (indicated in arcminutes on the con- 
tours in the figure). All of these contours, however, are 
defined for the same value of g'stCpix set by the governing 
parameters y and (Tb (or flfwhm)- In contrast to the opti- 
mal routine, a detection in the standard case is specified 
by not one, but two parameters - the pair {qst. Odd)- This 
embodies the fact that a detection must satisfy two cri- 
teria: a minimum flux density, ~ feCpix^dct/^pix' S-^d a 
minimum surface brightness, ~ Qstf^pix /fipix. In practice, 
the choice of values for qst and ^dct may be somewhat 
of a black art, but once made it specifies the survey's 
characteristic Af|J;j(z). For the ensuing examples, I make 
the choice motivated by the following considerations: Note 
that the signal-to-noise of a detection S/N = Qstx/Mnin- 
The detection criteria imposed as ^dct — V-^^min^pix may 
then be expressed in terms of the S/N: 

^ i (^) (^) ^''^ 

For reference, recall that in the optimal approach the pa- 
rameter Qopt was exactly the S/N. Now, at fixed signal-to- 
noise, a larger ^st leads to a smaller N^^ (i.e., 9det), which 
facilitates the detection of fainter point sources, because 
their flux density is buried in less noise (fewer pixels). On 
the other hand, a large value of qst disfavors finding low 
surface brightness objects; thus, a compromise is called 
for. One reasonable choice would be 6'dct ~ (l/2)(?fwhm, 
corresponding to a minimum detection S/N ^ 3, with 
Qst ~ 3 and 0pix/^fwhm ^ 1/2. I henceforth adopt these 
values for the following examples, which now completely 
specifies our detection routine. 

The dot-dashed lines in Figure 1 show the resulting 
standard mass detection curves. They all lie above the 
optimal detection curves, implying a lower overall sensi- 
tivity, as expected; and as in the previous cases, they fall 
with z. The low resolution examples in Figure Ic show 
a slight turn-down at low redshift, but under no circum- 
stances will they ever reach the origin at z = 0, as do the 
unresolved and optimal resolved detection curves: as al- 
ready mentioned, there always remains a non-zero detec- 
tion mass at low z in the standard case, due to the surface 
brightness cut. This constraint may be neatly summarized, 
using our earher result, as yo[Mf^^{z = 0),0] = y/27r (but 
it must be noted that this relies on our use of the self- 
similar cluster model) . The loss of close-by, low-mass ha- 
los is particularly noteworthy for the study of low-mass 
halos; to find them in an SZ survey, one of the important 
potentials of such efforts, will require special "tuning" of 
detection criteria, to more closely approach the optimal 
routine. 

5.2.2. Counts 

Our next goal is to use the detection mass to calculate the 
cluster counts and redshift distributions. It is worth not- 
ing in passing that one can envision several different kinds 
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of source counts: as a function of the value of ^dct, as a 
function of some aperture flux density (fixed or isophotal) 
or as a function of survey sensitivity, CTpix- Only the last 
one, however, is useful for optimizing an observing strat- 
egy, i.e., to answer the question of whether it is better to 
'go deep', or to 'go wide' when performing a survey. We 
are thus brought to consider in detail the detection mass 
as a function of detector sensitivity - M^*^(z; Cpix). The 
most direct efficient way to do this for a large number 
of different sensitivities is by returning to Eq.(p^), fixing 
Odct{— 1/26'fwhm) and then finding the detection mass as 
the root of the equation for each z. This avoids having to 
calculate the entire ^dot surface for each CTpix (y) just to 
extract a single contour. 

Using the result of this operation in Eq.(p^), we find 
the counts and redshift distributions displayed in Figures 2 
and 3 as the dot-dashed lines. Not surprisingly, the counts 
are lower at a given sensitivity than the corresponding 
optimal counts, and they are also slightly steeper. Similar 
remarks apply to the results of Figure 3. All of this is easily 
understood from the loss of low-surface brightness objects 
relative to the optimal routine. The essential conclusion 
concerning observation strategy remains the same: down 
to low flux densities, deeper integrations should yield more 
sources. 



6. Discussion 

The effects of resolving clusters must be properly modeled 
to understand the capabilities of possible ground-based 
surveys, as is clear from, for example. Figure 2: predicted 
counts are lower and steeper for resolved clusters relative 
to hypothetical SZ point sources. Besides lowering the ex- 
pectations for the number of detectable sources, these re- 
sults also suggest that deep integrations are more efficient 
than wide and shallow ones. The actual number of clus- 
ters expected for realistic ground-based performance are 
model dependent. For a self-similar cluster population, 
one could reasonably expect between 10-100 clusters/sq. 
deg. down to 0.1 mJy at 2 mm and with Sfwhm = 1 ar- 
cmin, as shown in Figure 2a. This number depends in 
addition on the source detection method employed: the 
standard routine counts may perhaps be considered re- 
alistic, while the optimal method counts indicate instead 
the best one could hope to achieve. At 1 cm, for equiva- 
lent sensitivity and at the lower resolution of 0fwhm = 2 
arcmin, one expects an order of magnitude lower surface 
density (Figure Ic). In sum, a square degree survey at 
2 mm could yield ~ 10—100 detections depending on 
the exact cluster model and the detection algorithm; a 10 
square degree survey at 1 cm to the same sensitivity (0.1 
mJy) could produce similar numbers. Both types of sur- 
vey may soon be achievable, with instruments similar to 
BOLOCAM (Glenn et al. 1998) or a detected interferom- 
eter array (Carlstrom et al. 1999) 



One of the primary interests of opening this new win- 
dow onto the Universe is to the search for high redshift 
clusters. The details of resolved cluster detection do not 
change the important and tell-tail difference between the 
redshift distributions in different cosmological models: the 
expected number of high redshift clusters is a sensitive 
function of JIq, as demonstrated by the redshift distri- 
butions given in Figure 3. Observations of such redshift 
distributions should prove a valuable tool for constraining 
r^o and for understanding evolution of the cluster environ- 
ment. 

There are several important issues that have not been 
dealt with in the present work. One concerns eventual 
source confusion, an effect that depends on the beam size 
and the exact value of the counts. This effect may very 
well be important even on arcminute scales, as noted by 
Aghanim et al. (1997). As these authors also point out, 
the issue is complicated by the fact that, due to the ex- 
tended nature of clusters, one must also contend with 
source blending. Detailed modeling of these effects really 
requires simulations. 

Another important issue not addressed in the present 
work concerns the question of radio source contamination. 
With sufficient frequency coverage, on can always identify 
SZ sources by their unique spectrum. Most often, though, 
spectral coverage is limited and contamination may be- 
come problematic. Its importance depends on the obser- 
vation frequency, and the counts at millimeter wavelengths 
are in fact a subject of current fundamental research; thus, 
the nature of contamination at in the millimeter is much 
more model dependent than in the centimeter. 

Finally, I note once again that the present work is 
based on a simple cluster model, because the principal 
motivation has been to understand the nature of resolved 
cluster detection by comparison to the more classic un- 
resolved case. Any attempt at a more exact examination 
of the number counts and redshift distributions requires 
more detailed cluster modeling. Such work would, in ad- 
dition, permit an interesting comparison of the relative 
efficiencies of SZ and X-ray observations to finding high 
redshift clusters, in practice. The SZ effect is clearly in- 
herently more efficient, but to really address this question, 
one should consider the actual achievable sensitivities of 
the two approaches. 

7. Conclusions 

There are clear and important differences in the con- 
clusions one draws concerning SZ surveys depending on 
whether clusters are considered as point sources or as ex- 
tended. For low resolution surveys, such as expected from 
the Planck Surveyor, most clusters will remain unresolved; 
however, when discussing the arcminute resolution more 
applicable to possible future ground-based surveys, we 
have seen that it is important to model the clusters as re- 
solved sources in order to properly understand the nature 
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of detectable objects. For a given sensitivity, high angu- 
lar resolution "resolves out" some clusters, lowering and 
steepening the final source counts. Relative to optimal re- 
solved detection, standard algorithms tend to in addition 
loose low mass, low redshift clusters due to their imposed 
surface brightness cut, further steepening and lowering the 
counts. With a fixed total observation time and a given 
frequency and angular resolution, we have seen that our 
results imply that deep integrations yield more objects 
than shallow ones covering a large area. 

Some important issues still to be explored concern 
the questions of source confusion and blending, and ra- 
dio source contamination. A detailed comparison of SZ 
and X-ray surveys would also be of interest, which im- 
plies more detailed cluster modeling than employed here. 

All the same, the numbers from the self-similar clus- 
ter model should be, within all the present uncertain- 
ties of these predictions, illustrative of what may be soon 
achieved from the ground. It appears that both in the 
millimeter and in the cm, ground-based SZ surveys could 
be capable of detecting up to ^ 100 clusters in total, a 
respectable statistical catalog. 
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